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The Quaternion Formulae for Quantification of Curves, 
Surfaces and Solids, and for Barycentres. 

By W. I. Stringham, Fellow of tlie Johns HopJdns University. 

In order to avoid a clumsy circumlocution, I have ventured to use the 
word quantification to denote in general that class of operations expressed in 
the several special cases by the terms, rectification, quadrature and cubature. 
Some of the quaternion formulae for quantification were given in a paper 
entitled " Investigations in Quaternions," communicated to the American 
Academy of Arts and Sciences., 9th January, 1878. They are here reproduced 
in a more general form together with the formulae for barycentres. 

Quantification. 

Let dM represent in general the element of arc, surface, or volume. Then 

if p = -^l/ {t) be the vector equation of any curve in space, where ^ is a scalar 

variable, 

dM-ds-Tfdt, (1) 

where p' = D^p =: the tangent to the curve at the point p. The formula for 

rectification, therefore, is 

szz^Tfdt. (2) 

The double area of the triangle, two of whose sides are p and p', is TVpp', and 
the element of area swept by p is 

dM =dA = ^ TVpp'. df . (3) 

Hence for quadrature of plane areas 

A = ^fTYpp'.dt. (4) 

This formula is sufficient to determine the area of a sector whose vertex is at 
the origin. In order to determine the area of any other sector (the surface 
being plane) whose vertex is S, it is only necessary to write 7t = p — S instead 
of p, 8 being the vector of the new origin with reference to the old. Thus, to 
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determine an area limited by a given chord, make h the vector of one ex- 
tremity of the chord and determine the limits of integration by assuming 
p =: 5 and p — ^ + >/, where yi is the chord in question. 

A surface in general is represented by the vector equation ^^z^it, u), 
where t^ u, are independent scalar variables. If at any instant u remain con- 
stant and t vary, p will describe a determinate arc upon the surface ; if t 
remain constant and u vary, p will describe another arc cutting the former at 

the point p. Hence 

p'l = D,p , p'2 = D,„p (5) 

will represent two tangents to the surface intersecting at their common point 

of contact. Therefore, 

dM—dS— T Vp'i p'2 . dtdu (6) 

will represent an elementally parallelogram upon the surface. Hence, for 

quadratures in general, 

8 = j'\lY^\^'^.dt,u, (7) 

where dt,u means dtdit, . The equation ^ z^ j^ [t , it) , for the special case of 
plane surfaces, may be written p = u^ {t) = un , the origin being on the sur- 
face. The formula last written then becomes 

S=^[u'-}ifTYrr'.dt, (8) 

where t' = D^r . This formula is a more general expression for (4) wherein 
the limits for u were and 1. It is evident that the choice of any other limits 
for u would determine an area lying between two similar curves. 

If, in the equation of a plane curve, p =: 4- (^) ? <|> be the angle of revolu- 
tion, and the curve be revolved about an axis a, (with the condition Ta = 1), 
then the element of area of the surface of revolution will be 

dM = dS = T Vap . d<p . Tp'. dt, (9) 

whence S = [<|>]* jTp'Vap . dt, (10) 

where p' = D,p as before. 

Let now u^rzp, where p:=xi^j ^) ^^^ z, t,u are independent scalar 
variables. When z varies between and 1, the extremity of q generates the 
solid enclosed by the surface whose vector generator is p. Since So'i^'jo's is 
the volume of the parallelopiped three of whose edges are g)\ , co'2 , o'g , then if 

DfCo = o'l , D„o = o'a , DjCO =: co'g , 
an elementary parallelopiped within the solid will be represented by 

dM=dV=S>o)\aW3-di,'i'',z. (11) 
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Hence, for cubatures in general, 

V=j'jj^a\J^a'^.dt,u,z, (12) 

and since o'l =: zp'i , (a'a = zp'2 , a'3 = p , 

.-. F=yra//Spp\p',.(??;,«, (13) 

where p'l = D^p , p'2 r=z D„p . Tliis formula enables us to determine the con- 
tents of sectors, or cones whose vertices are at the origin. For the cubature 
of a portion of the solid limited by a plane, instead of p write p — (^ , where h 
is the vector of some point in the given plane, and determine the limits of 
integration by making p satisfy the equation of the plane in question. 

The equation, p = ^^^ [t) , of a plane surface, with the condition that the 
surface be revolved about an axis a, will be sufficient to determine all points 
of a space enclosed by a surface of revolution. If Ta zr 1 and ^ be the angle 
of revolution, Vap is the projection of p on the normal to the plane of ap and 
Vap . d^ is identical with what ^'^du becomes in this case, that is, it is the tan- 
gent to a parallel of latitude. Then 

dMz^dV- Spp'Vap . zMtdzd^p , (14) 

p' = D;p . Hence, for solids of revolution, 

F=:-~[2V]J|j/Spp'Vap.fZ^. (15) 

The order of substitution of limits is 1 , z , afterwards 4) , . This may also 
be written 

V=\ lz^-]l I ^jTVpp' Vap . dt ; (1 6) 

for Spp'Vap = SVpp'Vap , and V (Vpp'Vap) =: ; therefore, 

S^Vpp'Vap = T^Vpp'Vap , or Spp'Vap = ± TVpp'Vap . 
In the formulae for cubatures, the limits and 1 for z give the contents of 
the space swept by the generating vector of the surface; any other limits 
would determine the contents of a space lying between two similar surfaces, 
i. e. of a shell. 

Barycentres. 
In Hamilton's equation* 

2m^ [e,. — e)-=zO, or 2???,.e,. = , (17) 

wherein the e and e,. are symbols of position in space, the s,. are vectors to the 
points e,. from e, and the m,. are scalar quantities, or weights, with which the 

* Elements of Quaternions, p. 89. 
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points «,. may be regarded as being loaded, — the point e is defined as the bary- 
centre of the system. 

Of the system represented by this equation the barycentre is evidently at 
the origin of vectors. A change of origin, referred to which y is the vector 
to the barycentre, gives _ 

2m, (e,, — y) =0, 
that is, 

7 — ^' — • (18) 

For a continuous, homogeneous mass, this equation assumes the form 

y- ^^.dM-^ ^clM, (19) 

where dM is the element of mass. It remains to substitute in this formula 
the values already obtained for dMva. the several cases above considered. 
For arcs dMrzz Tp'.t?^, (3), and 

y - JpTp'. dt -^ jTp'. dt. (20) 

For surfaces in general dM-=z TS^^i^<i.dt^u^ (6), and 

y = JJpTVp',p'2.fZif,« 4- JjTVp'.pV^!^,M. (21) 

For plane surfaces (8) this becomes 

y = y [^6«]i JtTVtt'. dt -r- \ lu% JtV^t'. dt, (22) 

^<T = p, t' rr D,Tr. For surfaces of revolution dM is (9), and 

y - JpTp'Vap . dt -^ jTp'Vap . dt, (23) 

a = axis of revolution, Ta zz 1. For solids in general dM is (11), and 

y= i i j (oQo\o'2Co's.dt,ti,z -^ j j (SG)\a'M'3.dt,u,z 

= \- MJ //pSpp'ip'2 -dt^u^^ \f-\l J/Spp'ip', , dt, u , (24) 

p'l ■= D(p, p'a = D„p. For solids of revolution dM is (14), and 

y- ^ [z^]J JpSpp'Vap.^?; - ^ [r]i, JSpp'Vap.(f^, (25) 

p' z=z D,p. This last formula [see (16)] may also be written 

r = \ M' /pTVpp'Vap .dt-^^ \/-]l jTVpp'Vap . dt. (26) 

It is to be noticed that the formulae (24), (25), (26) give the barycentres 
of the entire solids, or of shells, according as the limits assumed for z are 
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and 1, or some other. Suppose the limits to be z and z + e, wliere e is infini- 
tesimal. If the terms containing the second and higher powers of s be 
neglected, there will be in the numerator the factor (z + ^Y — z^=:4z^£, and 
in the denominator, (z + e)^ — z^zz 3z^e ; and (24) will reduce to 

y^z i j pSpp'ip'a.f?;^,^ -r- j j Spp'lj^'2.f?^5,^f, (27) 

the barycentric vector to an intinitely thin shell. The relative thickness of the 
shell, as shown by the factor z in the equation, o = z^ {t,ic), of the shell, will 
be determined by the position of the origin of the generating vector. 

AjfjMcations. 

The ellipsoid will afford a convenient illustration of the application of the 
above methods. Its equation is 

p-=a cos X -{-a sin X , 
where cr = /3 cos y + y sin y , 

and a, (3, y are the three principal semiaxes. By differentiation and reduction 

TVp'ip'2 ^ as sin x Vl — sin^^o cos'a; , 

b^o^ be 

where ar^ Ta , h-zzTj3 , c=: Ty , s = Ta', and sin^io rz 1 ^~^, or cos it; = — . 

CI s cts 

Hence by (7) 

Szrz] I sin X Vl — sin^ w cos'^ x . dx , y . 

J cos w J "^ 

Let sin ■?; = sin 10 cos a; , cost) — Vl — sin^'w cos'^a; . Then 

„ r 2bo r 2 , fbe \_2v + sin 2t)]S , , 

When the surface is one of revolution, say the prolate ellipsoid, then h^c=zs, 

C^ — IP' 

sin^w rz — -^ , and the expression (28) reduces to 

,Sz=^-;^J^(2« + sin2o)]SU. 

The order of substitution of limits is v , , afterwards i/ , . 

By easy transformations it will be found that Vpp'i = Vaa , Yaa' = Yf^^y, 

and therefore 

Spp'ip'2 = abc sin X . 

Hence, for the ellipsoidal shell [see (13)] , 

V=^ [z^]i jysin x.dx,^ = "g" [z^t/ cos a']; \'i\l. 
The order of substitution of limits isa;,0; j/,0; l,z. 
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The value of pSpp'ip'2 is ahc {a sin a; C0Sor + a sin^^), and its complete in- 
tegral is 

"- [2ay sin^ x + (/? sin if — y cos ?/) {2x — sin 2x)'}l \ K . 

Hence the barycentric vector of the ellipsoidal shell is 

— 3 [s* {2ai/ sin^x + (/9 sin y— j- cos y) {2x — sin 2.-c)}]J | ? | J 

^ ~ lQ[fymsx][\i\l ' 

The order of substitution of limits is important ; it is as above in the expres- 
sion for V. If the shell be infinitely thin, in this case an ellipsoidal Ohasle- 
sian shell, the application of formula (27) gives 

— 2 [2ay sin^ x + (;9 sin y — yco» y) {2x -~ sin 2a;)]J 1 1 

^~ 4[2/cosa;]5lS 

For the ellipsoidal solid y is the expression last written multiplied 

3 
by -J- . The barycentric vector to the half solid, bounded by the plane 

3 
of a(3, will easily be found to have the value -r- y, and that of one-eighth 

of the solid, bounded by the planes a/3, f3y, ya, to have the value -k- (a + Z^ + y). 



It is worth Avhile to observe, in this place, that from the equation 
p=: q — to , (where p is the generating vector to a family of surfaces, and w is 
a scalar parameter), which Tait* has assigned as the equation of a volume, 
would be deduced substantially the same results as the foregoing. Let T^ 
represent the length of the radius vector to the surface which bounds the 
solid, and let q be written as a function of t , u , 6 ; then if t,u be taken as 
the variables which define the surface, T^' and UV^' will be functions of t, u 
exclusively, and will turn out to be the angle of q. Thus the equation of 

the solid may be written 

q = ^{t,u,e), 
or 

Yq = X (^) ^0 ^i'l ^ • 
This vector equation will be satisfied for every point within the solid ; and it 
has the same form as the equation p = 3% {t, u), made use of in the foregoing 
discussion. 

^Treatise on Quaternions, p. 61. 
April 15th, 1S79. 



